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Abstract 

Following the anomalous like-sign dimuon charge asymmetry measured by the 
DO collaboration at the Tevatron collider we discuss the implications of large CP 
violation in B^ s mixing for Supersymmetric (SUSY) Standard Models, focussing 
on those models which allow a family symmetry and unification. For the Minimal 
Supersymmetric Standard Model (MSSM) we show that it is only possible to 
account for B s mixing and CP violation at the expense of large squark mixing 
which would require a new approach to family symmetry models. In order to 
describe both B s and B^ mixing and CP violation we are led to consider SUSY 
models with Higgs fields transforming as triplets under a family symmetry. We 
describe a realistic such model based on A27 family symmetry in which tree-level 
exchange of the second Higgs family predicts B s and Bj mixing and CP violation 
in good agreement with a recent global fit, while naturally suppressing flavour and 
CP violation involving the first and second quark and lepton families. 
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1 Introduction 



The Standard Model (SM) has provided a remarkably successful description of quarks 
and leptons from its inception in the 1960s until the end of the last millennium. In 1998 
the discovery of neutrino mass and mixing demanded new physics beyond the SM for 
its explanation (for a review see e.g. p]). The picture which has emerged in the lepton 
sector is consistent with three neutrino mass and mixing described by a PMNS matrix U, 
although its origin remains unclear and U e ^ is so far unmeasured. By contrast, despite 
intense experimental and theoretical scrutiny, there has been no firm evidence of any 
new physics in the quark sector, with the CKM picture of CP violation, summarised 
by the Unitarity Triangle, becoming ever more precisely determined [2]. Yet despite 
this progress, some cracks have begun to appear in B physics which may call for new 
physics to describe CP violation beyond the CKM matrix. It is worth recalling that 
CP violation is predicted by the SM to be very small in B — B mixing, well below the 
Tevatron sensitivity, due in part to the small phases of the relevant CKM elements. On 
the other hand new physics can compete with the SM box diagrams, in principle with 
large new CP violating phases, giving much larger CP asymmetries in B — B mixing 
than predicted by the SM, rendering it observable at the Tevatron. 

Recently the DO Collaboration has reported evidence for CP violation in the like-sign 
dimuon charge asymmetry [3] 

N ++ — N — 

At = — ^rx — = -(0.957 ± 0.251 ± 0.146) x 10~ 2 , (1) 

where (N b ) is the number of events with b (b) containing hadrons decaying 

semileptonically into fi + X (yrX) . The DO result is 3.2a away from the standard 
model (SM) prediction (—2.3 ±0.5) x 10 -4 jl]. The CDF [5] measurement of A b h using 
only 1.6 fb _1 of data, as compared to the DO 6.1 fb _1 data set, has a central value 
which is positive, A b sl = (8.0 ± 9.0 ± 6.8) x 10 -3 , but is still compatible with the DO 
measurement at the 1.5a level because its uncertainties are 4 times larger than those of 
DO. Combining the the DO and CDF results for A b sl , one finds A b sl ~ -(0.85±0.28) x 10" 2 
which is still 3cr away from the SM value. 

The interpretation of the observed CP asymmetry is in terms of the production of 
BB meson pairs followed by their subsequent oscillation and semi-leptonic decay where 
the charge of the final state lepton effectively tags whether it is a b or b quark which 
decays. Thus the dilepton asymmetry can be written as 



Ab ^ N(BB) - N(BB) = p^=-p^T b 
sl N(BB) + N(BB) p b^b + p b^b " 

p b^b Pb^b 

The measured asymmetry at the Tevatron is interpreted as a linear combination of the 
asymmetries a s f in Bd and B s oscillation and decays [3], 

A b sl = (0.506 ± 0.043)4 + (0.494 ± 0.043)4, (3) 
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where the "wrong charge" asymmetries are 
q = T{B q ^^X)-T{B q ^^X) 

flsl r(B q ^ii+x)+r(B q ^^-x) • 1 ' 

The current experimental values of the separate asymmetries are afj = —(0.47 ±0.46) x 
10~ 2 [6] and a s sl = -(0.17 ± 0.91 ± 0.15) x 10 -2 [7] which, while being consistent with a 
negative of order one per cent, are also consistent with zero as well, and so do not 
shed much light on which of the two separate asymmetries is responsible, however there 
is apparently a mild tendency for both of these asymmetries to be acting together. 

In addition to these measurements, DO have reconstructed B s — > J /ip<t> decays and 
have measured the time dependent asymmetry parameter = — sin0£ s , where <pB s 
is the phase of the B s — B s mixing matrix element M\ 2 = |M s 12 |e^ Ss , and finds a 
discrepancy with the SM prediction of ~ at the level of 2.1a [8]. A recent 
preliminary CDF analysis based on 5.2 fb _1 of data, finds which is consistent with 
zero but has a central value of ~ 0.5 [H] . The predictions for in various models 
and the correlations between this and other observables as a means of discriminating 
between these models has been comprehensively studied [10] . Neglecting the small SM 
contribution to S^, the following model independent relation holds between and 
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■SV>0 (5) 



|Mi 2 , 

Following the recent DO results, it has been shown that it is possible to fit and 
fr° m t ne latest data by assuming that new physics contributes significantly to the 
mixing matrix element Mf 2 and also by allowing the decay matrix element r* 2 to float 
[12] . More precisely the authors in [12] perform a global fit of all experimental mea- 
surements, including the recent DO asymmetry results and the recent CDF preliminary 
results for B s — > J/ip4>, allowing the two SM decay matrix elements Y\ 2 to float, while 
allowing for new physics to contribute to both mixing matrix elements Mf 2 which can 
be parameterised as: 

Mf = Mf 5M (l + h q e w "). (6) 

Using the convention where = Im(— r^ 2 /M^ 2 ) with the dominant real parts of 
and M( 2 being positive in the SM, gives, 



l sl ~ 



|r* 2 | h s sin6 s 
Ml 2 ' SM \l + 2h s cos9 s + h 2 '' 



(7) 



where we have neglected the small phase in the SM matrix element, S ~ —0.01. How- 
ever the corresponding SM matrix element in the B^ sector is, M^ 2,SM = |Mj 2 ' 5M |e 2l/3 , 
with j3 ~ 0.38, which will contribute significantly to the phase of M\ 2 . 
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The global fit [12] includes the measured B s and Bd mass differences AM S and AM d , 
the measured time dependent asymmetries S^, S^k (which determines the unitarity 
triangle angle /3), the CP asymmetries a sl ' s as well as the CKM parameters p, fj, while 
allowing T\ 2 to vary from its SM predicted value. The best fit value for is about 
twice as large as the SM prediction [3], with similar results obtained in [T3] where 
the implications of such a large value are discussed. The decay matrix element T* 2 is 
proportional to the square of a tree-level SM amplitude proportional to V c b arising from 
W exchange so it is challenging to understand why the best fit value should be so large. 

The best fit points for new physics contributions to the matrix elements are [T2] 



Note that there are two different best fit points for (h s ,6 s ), but only a single best fit 
point for (hd, Bd) which is only preferred from a zero value at a confidence level of order 
1.5cr. All points have the angle 9 q in the third quadrant where both sin 6^ and cos8 q are 
negative, resulting in enhanced negative values of a q sX from Eq|7j Although the precise 
best fit points must be regarded as indicative values with rather large error bars, the 
fit is quite robust with h s = hd = disfavoured at 3.3cr [12J. There have already been 
several attempts to explain the recent data [HI |15l [161 H3 HH UHl ED] (see also [21J). 

In this paper we consider the implications of large CP violation in B mixing for SUSY 
Standard Models (for a review see e.g. [22]) focussing on those models which include 
a family symmetry and allow for unification. We discuss two distinct possibilities. For 
the Minimal Supersymmetric Standard Model (MSSM), discussed in Section 2, we show 
that it is only possible to account for B s mixing and CP violation at the expense of 
large squark mixing which would require a new approach to family symmetry models. 
Moreover this approach cannot account for Bd mixing and CP violation. In order to 
describe both B s and Bd mixing and CP violation, in Section 3 we are led to consider 
SUSY models with Higgs fields transforming as triplets under a family symmetry, where 
tree-level exchange of the second Higgs family may readily account for all the data. We 
describe a realistic such model based on A27 family symmetry [23] which predicts B s and 
Bd mixing and CP violation in good agreement with the best fit point I, and naturally 
leads to small effects in K° mixing and other flavour violating processes. 



2 Large CP Violation for B s mixing in the MSSM 



As discussed in [23] the MSSM contributions to CP violation in B mixing arise domi- 
nantly from box diagrams involving down squarks and gluinos. The B q mixing matrix 
element can be written as the sum of the SM box diagrams and the SUSY box diagrams, 



(h d , 9 d ) ~ (0.25, -tt), (h s , 9 s ) f ~ (0.6, — tt), (h s , 9 s ) n 




(8) 




12, SUSY 



(9) 
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where from Eq(6]we identify 

j^l2,SUSY 

h e i9q = - CIO) 

where q — s,d. Using the matrix elements in [21] with the updated parameters in [17] 
we find, in the mass insertion approximation (see [25J and references therein), 

h*»' - f^^^^Bft^ + ^^-lMK^M^^ (11) 
V m q J 

hde m +m w f 500^eV\ 2 [n4 + _ 446Q {{ s$ i)ll{§ $ i)rr)] , 

where (<^) LL = (V Dz ,m|Vj i ) 3i /m| and (S^rr = (U Dfl m|Vjj 3i /m?, with m 9 ~ being a 
typical squark mass (assumed to be degenerate with the gluino g) where mj~ is the left- 
handed (L) squark doublet mass squared matrix, m 2 ^ is the right-handed (R) down-type 
squark mass squared matrix, and Vd l and Vd r are the unitary matrices that diagonalise 
the down-type quark mass matrix M d , namely VD L M d V^ R = diag(md, m s , m;,). This is 
summarised by the statement that the LL and RR mass mixing between down squarks 
of different generations in the super CKM basis is the source of the flavour and CP 
violation [22]- We have not included the contributions from LR mass mixing which are 
tightly constrained by b — > sj. 

Comparing the SUSY predictions in EqjTT] to the best fit points in Eq|HJ it is clear 
that the values of h s ~ 1 could either be achieved, assuming squark and gluino masses 
of about 500 GeV, by (<5f 2 ) LL ~ (5l 2 ) RR ~ 0.05 - 0.1 or (<&) i£ < (5l 2 ) RR ~ 0.3 - 0.6 
or (fi^RR ^ ($32) ll ~ 0.3 — 0.6. These represent quite sizeable squark mixing angles, 
which run into conflict with grand unified theories (GUTs) based on gravity mediated 
SUSY breaking. To see this, it is worth bearing in mind that renormalisation group 
(RG) running from the high energy GUT or Planck scale to low energies tends to 
increase the diagonal squark masses rriq by about a factor of 5, while not enhancing 
the off-diagonal squark masses [22] , so the high energy (5^ 2 )ll,rr parameters need to 
be 25 times larger than these low energy values which is not possible (they can at most 
only be of order unity). For some gauge mediated SUSY breaking scenario, where the 
messenger scale is below the GUT scale, the effect of running is reduced so it may be 
possible to achieve these low energy values. Another constraint is that, in the framework 
of GUTs, there is the danger of running into conflict with the bound on (<5| 2 )ll < 0.12 
from r — 7- /i7 [25J since the slepton masses are only enhanced by about a factor of 
2 in running from the GUT scale to low energies. In the context of SU(5) GUTs 
the low energy (S^ 2 )ll parameters are constrained by the high energy requirement that 
($32) rr T = ($32) ll T < 0-48 [25], which implies the low energy constraint (6$ 2 )rr < 0.02. 
From the point of view of family SU(3) symmetry models (see [26J and references therein) 
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small squark mixing parameters are also expected {S 32 )rr ~ 10~ 3 . It seems that the 
data is not favoured by conventional SUSY GUTs and family symmetry. 

Suppose we abandon all pre-conceived prejudices about SUSY breaking, GUTs and 
family symmetry, but continue to assume that the MSSM is the only source of new 
physics. Then we can ask if it is possible for the MSSM to describe the observations 
and if so then what the data is telling us about squark mixing. For the reasons outlined 
above, from the point of view of the MSSM, it is desirable for the squark mixing to be 
as small as possible. Therefore we shall consider the smallest mixing describing the data 
given by solution I with (S^ 2 )ll ~ (^32) RR ~ 0.05. Taking into account the RG running 
up to the (unknown) SUSY breaking messenger scale, this still suggests a high energy 
theory capable of giving quite large (2,3) mixing in the squark sector. Following this 
reasoning we are led to consider high energy Yukawa matrices in the quark sector which 
have a democratic structure in the (2,3) sectors, 

Y&A ~ Y(U ~ (J })■ (12) 

This ensures that even approximately diagonal squark mass squared matrices would 
generate large squark mixing in the super CKM basis. However the small CKM angle 
#23 ~ I #23 — 6*23 1 "would then require an accurate cancellation. To enforce this (ap- 
proximate) cancellation in a natural way we shall require Y% 3 ) oc Y& ^ and rank one 
sub-matrices to achieve the (2,3) quark mass hierarchies. Both these features could be 
achieved by an 577(3) family symmetry under which the left and right-handed quarks 
transform as triplets Qi, U?, D\ ~ 3 and which is spontaneously broken by an anti-triplet 
flavon with aligned vacuum expectation values (023) ~ (0, 1, 1) V, where phases have 
been suppressed (see [26]). Then the (2,3) block could be generated from leading order 
operators of the form, dropping coefficients, 

H u Q l( pi 3 U^ 23 + H d Q^D^ (13) 

which, after the two Higgs doublets of the MSSM H u and H d acquire their VEVs, 
implies #23 ~ and m s <C m^, m c <C m t , at leading order. This differs from the usual 
SU(3) models [2S] by the absence of the flavon 3 with VEV (fa) ~ (0,0, 1)V 23 . This 
democratic (2,3) structure could be plausibly be extended to the charged lepton sector 
as well, 

*? 2 ,3)~(j I), (14) 
resulting from the leading order operator, 

H d L i( f>i 3 E°4 3 , (15) 

which implies maximal charged lepton mixing in the (2,3) sector. In order to achieve 
maximal (2,3) physical lepton mixing we require the light effective neutrino Majorana 
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mass matrix to be approximately diagonal. For example this could be due to a type 
II see-saw mechanism via a sextet flavon A 4 - 7 with an approximately diagonal VEV 
(A*- 7 ) ~ diag(0,a,6) arising from operators of the form |27j . 

H u H u LiLjA ij . (16) 

This then leads to maximal atmospheric neutrino mixing coming from the charged lepton 
sector. 

In the above approach we are essentially saturating the high energy limits with 
(<&) lljir ~ (d^LLJiR ~ 1 leading to some apparent tension with the previously men- 
tioned r — >■ fi'y limit (5| 2 ) ^ T < 0.48. However there is a very large compensating effect 
from RG running due to the large charged lepton Yukawa couplings which will tend to 
reduce the magnitude of the slepton masses, including the off-diagonal ones. To leading 
log, the correction will be 

A(5 3 e 2 )LL~-^ln%^~-0.4 (17) 

87T 2 M W 

effectively relaxing the tension and allowing (S^ll ~ 1- Nevertheless, r — > /ry might 
be expected to be not far below its current limit. Such Yukawa induced RG corrections 
will also be present for the squark sector, but there the dominant suppression is coming 
from QCD enhancement of the diagonal squark masses. 

Turning to the less statistically significant evidence for Bd mixing, comparing the 
SUSY predictions in Eq.ll to the best fit points in EqjSJ it is clear that the values of 



h d ~ 0.25 could either be achieved by (5 31 )ll ~ ~ 0.8 x 10~ 2 or (S^n < 

{$3i) rr ~ 5 x 10 -2 or (S^rr <C (^fi)ii ~ 5 x 10~ 2 . From the point of view of 
conventional family SU(3) symmetry models the above squark mixings are much larger 
than the expected values (S^rr ~ 10~ 4 [26]. Although these low energy squark mixings 
look more modest, they must originate from high scale squark mixings which are many 
times larger than these values (25 times if the high scale is the GUT scale). Again 
we consider the smallest mixing case corresponding to (S 31 )ll ~ (831) rr ~ 0.8 x 10~ 2 
corresponding to high scale values perhaps of order A ~ 0.2. This suggests a model 
with (1,3) quark mixing angles of order A in the up and down sectors, 8™ 3 ~ 8f 3 ~ A as 
compared to #13 ~ A 3 , which again demands some natural cancellation mechanism. In 
order to achieve this we may extend the SU(3) approach above, by introducing using 
the anti-triplet flavon cf>\ 23 with aligned VEV (0123) ~ (1, 1, 1)^123, where phases have 
been suppressed (see [2S]), which introduces the additional operators 

H U Q^\ 23 U^ 23 + H u Q l( p 23 U^\ 23 + H d Q l( p\ 23 D^ 23 + E A Q^ 23 D)<\? X23 . (18) 

The combined effect of the operators in Eqs. T3fT8 is to yield the quark Yukawa couplings, 



assuming V123/V23 ~ A, ignoring phases, 

/OA A \ 

Y u ~Y d ~ A 1 + 2A 1 + 2A . (19) 
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After the maximal angle (2,3) rotations, the remaining matrices are diagonalised by 
separate and equal (1,3) rotations of order A in both the up and down sectors which 
cancel as required. However this results in first family eigenvalues of order A 2 in both 
the up and down sectors leading to the relations m u /m t = m^jm^ ~ A 2 which are 
in strong disagreement with their observed values. This result could in principle be 
evaded when other operators or corrections are included as required to describe the 
second family quark masses m c , m s , but apparently only at the expense of an unnatural 
cancellation of the leading order first family result above. Therefore we are not able to 
accommodate large and natural (1,3) mixing of order A in the up and down sectors as 
suggested by the Bd mixing fits. In general we do not see how to achieve this in a natural 
way in the MSSM including family symmetry and GUTs. We conclude that either the 
weaker requirement of non-standard Bd mixing be discarded, which as remarked in the 
Introduction is only required at a confidence level of order 1.5a, or one must abandon 
the MSSM as a natural explanation for CP violation in both B s and Bd mixing within 
the framework of family symmetries. 

3 Large CP Violation for B s and Bd mixing with 
three families of SUSY Higgs 

It is obviously easier to account for large CP violation in B s and Bd mixing by the 
tree-level exchange of some new heavy particles than if they only appear in a one-loop 
diagram, as is the case for squarks and gluinos in the MSSM which only appear in the 
SUSY box diagram considered in the previous section. Indeed many of the approaches 
suggested in [TH [151 HSl EG] are based the tree-level exchanges of some new particles. 

A good example of this approach is the suggestion [15J that CP violation in B s mixing 
may be due to the tree-level exchange of a new neutral spin- boson H° A = (H°+iA°)/V2 
which is the neutral component of a Higgs doublet where H° does not get a VEV at 
leading order so that H° and A have approximately the same mass Mh = Ma- The 
Yukawa couplings of H® to b, s quarks in the mass eigenstate basis are given [15] , 



Tree-level H® exchange gives rise to the following operator which contributes to B s 
mixing, 



ivi H 

One may readily extend this idea to allow for couplings of H® to b, d quarks in the 
mass eigenstate basis, 



H d (y bs b R s L + y sb s R b L ) + H.c. 



(20) 



Ml 



(b R s L )(b L s R ). 



(21) 



Hl(y bd b R d L + y d bd R b L ) + H.c. 



(22) 
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Tree-level H^ exchange then gives rise to the following operator which contributes to Bd 
mixing, 

y -j^(b R d L )(b L d R ). (23) 

Following the approach of the previous section, the B q mixing matrix element can be 
written as the sum of the SM box diagrams and the Higgs tree-level exchange diagrams, 



M 12 = M 12,SM + M 12,Hig gS (24) 

where from Eq(6]we identify 

hq e " = ^12,SM ' (25) 

where q — s,d. Using the matrix elements in [17J we find, 

ft^«« » -3.4 x loWi ( 2 °^ eY ) 2 ■ (26) 



Comparing the Higgs predictions in Eq,26 to the best fit points in EqjSJ it is clear that 



the values of h s ~ 1 could be achieved, assuming Higgs masses of about 200 GeV, 
by Ubs ~ Usb ~ 10~ 3 , depending on the phases. Similarly, hd ~ 0.25 could be achieved, 
assuming Higgs masses of about 200 GeV, by ~ yab ~ 10 -4 , depending on the phases. 
In both cases the required Yukawa couplings are proportional to the Higgs mass. 

Of course such tree-level exchanges must be kept under control so that they don't 
induce too much flavour changing in other places where the constraints are more severe, 
especially involving the first two quark and lepton families as is the case for example with 
ex- To overcome these challenges, it has been suggested that the hypothesis of Minimal 
Flavour Violation (MFV) could be extended to the two (or more) Higgs doublet model 
[T7j . Here we shall follow a different approach, namely to use the idea of family symmetry 
to control the magnitude of the flavour changing Yukawa couplings in the framework of 
a model with three families of SUSY Higgs doublets. 

The basic idea we shall discuss is very simple, namely that there are three SUSY 
Higgs families which form triplets under some family symmetry group, just like the 
three families of quarks and leptons. The three SUSY Higgs families can be written 
as i?", Hf where the index i = 1, 2, 3 labels the three copies of the two MSSM Higgs 
doublets H u ,H d . The idea is that the three families of quarks Qi^Uf^DI and Higgs 
if", Hf all transform as triplets ~ 3 under some family symmetry group. Yukawa 
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couplings are forbidden in the exact family symmetry limit but are allowed when the 
family symmetry is broken by flavon VEVs. The idea is that the family symmetry 
leads to a doubly hierarchical structure as follows. The third family of Higgs couples 
more strongly to quarks and leptons than the first and second Higgs families. The 
third family of Higgs also couples more strongly to the third family of quarks than the 
first and second families of quarks. This leads to the observed pattern of quark masses 
and mixings due to the third Higgs family couplings, with suppressed flavour changing 
couplings from the first and second Higgs families. 

A model with three families of SUSY Higgs which controls the flavour changing using 
a family symmetry has recently been discussed [23]. The model [23] is based on a A27 
family symmetry and leads to a successful description of all quark and lepton masses 
and mixing, including tri-bimaximal lepton mixing [28J. However only the dominant 
third Higgs family couplings were considered [23J and flavour and CP violation arising 
from the subdominant first and second family Higgs couplings were not considered. In 
the following we shall revisit this model, focussing on the quark sector, assuming exactly 
the same particle content and symmetries as in [23] , but including the effects of different 
field orderings and contractions not previously considered. 

In the considered model [23] the quarks and Higgs transform as triplets under a 
A 2 7 family symmetry, which is broken by anti-triplet flavons 03,03,023,0123 which de- 
velop aligned VEVs, dropping phases in the following discussion for simplicity (they are 
recovered in Appendix |Al, 



(0 3 ),(0^)oc ( 1 ) T , (0 23 )oc ( 1 1 ) T , (0 123 )oc ( 1 1 1) T . (27) 

The leading order down-type quark Yukawa couplings result from the following flavon 
couplings, suppressing coupling constants, 

j^- [ (g03)lp C 03)l(^)l + (g023)l( J D C 023)l(^)l 

+ (Q012 3 )lp C 023)l(^03)l + (Q023)lp C 0123)l(^03)l ], (28) 

with similar couplings generating the up-type and lepton Yukawa couplings. The vacuum 
alignments in Eq]27]then imply that only the third Higgs family couples to quarks and 
leptons with the leading order Yukawa matrix, defined by Y^ 3 QiDjH^ } given as, 

/ 4 e\ \ 

yf ~ \ 4 4 4} (29) 

\4 4 i / 

where the expansion parameter e d ~ 0.15, where we have assumed [23] Q 

«>„<*>.., %> we , (30) 



M h M d M d } l,i 

iMore precisely « ^ « 0.8 [23] 
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A similar Yukawa matrix arises for the up-type quarks but with a smaller expansion 
parameter e u ~ 0.05. When supplemented by additional corrections [23], such Yukawa 
matrices, after the third Higgs family develop their VEVs, have been shown to provide 
a successful description of quark masses and mixing [29J . 



The couplings in Eq 28 are controlled by the additional symmetries of [23J and are 
generated by the exchange of heavy messenger particles of mass and Mh, as indicated 
in Figjlj The quark singlet messengers with masses are assumed to be lighter than 
the quark doublet messengers of mass Mq, and hence give the dominant contribution 
with M u w 3Md being responsible for ~ 3e u . The assumed messenger sector also 
ensures that 0§ couples directly to the Higgs doublets with the group theory contractions 
giving the A27 singlet as indicated by the subscript in (H d (f)^)i. This implies that only 
the third Higgs family couples to quarks. 

In Figj2] we display other messenger diagrams, not considered in [23], allowed by 
all the symmetries of the model and involving the same flavon fields but contracted 
differently. These operators will contribute at a suppressed level due to the (assumed) 
heavier primed messenger masses. The diagram on the left in FigfJ] will lead to the 
following additional operators, suppressed by two primed messenger masses, 

3^(^3)1(^3)1 + (Q03)l(^23)l(^ 2 3)l 



M d M> d M> h 



+ (d^)!^)! + (Q^i^ua)^^)! ], (31) 

with similar couplings for up-type quarks and leptons. The diagram on the right in Fig|2] 
will lead to the following additional operators, suppressed by three primed messenger 
masses, 

— \— [ (Q03)lp C ^)l(^3)l + (Q023)lp C ^)l(^23)l 
M d M h 

+ Wl23)lp C 03)l(^23)l + (Q023)lp C 03)l(^123)l ], (32) 

with similar couplings for up-type quarks and leptons. Other messenger diagrams with 
two flavons along the Higgs line will not introduce new operator structures but will 
only change the overall coefficient of the operators in Eq j3~Tf32"| with one power of M' d 
being replaced by M' h . Additional operator structures are also present in which the 
singlet contractions (represented by the subscripts) are replaced by the one dimensional 
representations 1', 1" present in A27, as discussed in Appendix |Aj However we assume 
that messengers in other one dimensional or higher dimensional representations are 
absent or very heavy. 

The new couplings in Eqs, 31~j32"| together with the vacuum alignments in Eqj27| allow 



the first and second Higgs families to couple to quarks and leptons with the leading order 
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H 

Q 



M d M d 



Figure 1: The messenger diagram for the down- type quark sector as assumed in [53j where <fr represents 
any of the flavons ^3, (^23, 0123 allowed by the symmetry. Analogous diagrams are present in the up- type 
quark and lepton sectors. 



H 1 



Q 



D' 



M' d M d 



Q 



Ml 



D' 



M' r , Mi 



Figure 2: Other possible messenger diagrams involving the same flavons as in FigJlJ Additional dia- 
grams with two flavons along the Higgs line are also possible but we assume that diagrams with three 
flavons along the same line are suppressed due to higher dimensional messenger representations being 
absent or very heavy. The primed messenger masses are assumed to larger than the unprimed ones, 
leading to the diagram on the left being suppressed compared to Fig[l] and the diagram on the right 
being even more suppressed. 
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Yukawa matrices, defined by Y^ k QiDjH^, given as, 














Y. 




4 a d 



(33) 



where we have used the expansion parameters in Eqj30} together with the suppression 
factors, 



ad 



Md 
M d 



a h 



Mh 

K 



Eq 33 shows that the second family Higgs H 2 



(34) 



couples more strongly to quarks than 
the first family Higgs Hf whose effects can be approximately ignored. Assuming as a 
leading order approximation that only the third Higgs family develops VEVs, which 
is reasonable since the VEVs are radiatively generated as a result of the large third 
family Yukawa coupling in Eqj29] of order unity, then we may diagonalise the quark 
mass matrix resulting from Eq.29 by small angle rotations 6^ ~ e 2 ,, 9f 3 ~ e\, 0f 2 ~ e^, 



to go to the diagonal down-type quark mass basis d,s,b. The neutral component of 
the second family Higgs H 2 ° then has a Yukawa coupling matrix (the first equation in 



Eq.33) which, when rotated to the down quark mass basis, takes the leading order form, 



c& d 

Y dl ss ~ I 4 4 4 a d I a d a h , (35) 

-3 ,2 

d fc cZ 




in the convention where the rows correspond to dL,si,bL and the columns correspond 



35 



arc 



to d}^, sr, bji- Clearly the flavour violating couplings of H^ involving bn in Eq 
relatively suppressed by a factor of compared to those involving as can be under- 
stood from Fig|2j Note that there are no cancellations of the couplings in Eq 35 in the 



d,s,b basis since each operator in Eqs 28pTp2 has an independent order unity coeffi- 



cient which has been suppressed for clarity The flavour violating couplings in Eqs. 20p2 



can then be read off from Eq, 35 , where we identify H® = H 2 



Vbs ~ 4 a d a h, Usb ~ 4 a d^h, Vbd ~ 4 a d a h, Udb ~ 4 a dOt h - (36) 



With the phases included, as discussed in Appendix [AJ the couplings in Eq|36] lead 
to CP violation with, 

\ ybdV *>\ . £ ^2x lO" 2 , (37) 



\ybsy; b \ 



which leads to the prediction, using Eq.26 

hd 1 
K ~ 2' 



(38) 
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in good agreement with the best fit point I in Eq|8j This prediction is independent of 
the messenger masses and hence is equally valid for other messenger diagrams with two 
flavons along the Higgs line. 

Assuming best fit point I the required Higgs mass may be related to a d , ah using 



Eq.26 



Af, 



H \ ~in3„.3„2 



200 GeV 



10^X- (39) 



For example ad ~ ah ~ 1/3 would require Mh ~ 400 GeV. This choice of couplings 
would then imply, 



y sb ~ e 2 d a d a h ~ 2.5 x 10~ 3 , y db ~ e d a d a h ~ 3.7 x 10~ 4 . (40) 



The hierarchical structure of flavour changing couplings in Eq,35 suppresses the 
contribution of second family Higgs exchange to K° mixing, which is described by a 
very small Yukawa coupling 

Vsd ~ Vds ~ 4a d a h ~ 8 x 10~ 6 , (41) 

where we have assumed a d ~ Qh ~ 1/3. The charged leptons are expected to have 



flavour violating couplings with a similar structure to Eq. 35 , with the hierarchical struc- 



ture again leading to suppressed lepton flavour violating processes with 

2/ M e ~ Hep ~ Vsd ~ Vds ~ 8 x 10" 6 . (42) 

The diagonal coupling of the second family Higgs to muons is also quite suppressed, 

Vim ~ e d a d a h ~ 6 x 10~ 5 , (43) 

leading to a negligible contribution to Br(B s — > /i + /i~). For example, assuming td ~ 
0.15, ad ~ ah ~ 1/3 and Mjy ~ 400 GeV using the results in [15] we estimate that the 
contribution from second family Higgs exchange to the branching ratio is ABr(B s — > 
// + /x") ~ 3 x 10~ 13 which is negligible compared to the SM prediction of about 3.7 x 10~ 9 . 
This contrasts with other non-standard Higgs models [T5"j ITT] which tend to predict 
Br(B s — > yU + /i~) larger than the SM value and close to the current experimental limit 
of about 5.8 x 10~ 8 . 

Similarly the diagonal coupling of the second family Higgs to taus is, 

Vrr ~ e 2 d a d a h ~ 2.5 x 10~ 3 , (44) 

leading to a new contribution, assuming the same parameters as above, ABr(B s — > 
t + t~) « 4.5 x 10 -10 , which is also somewhat below the SM prediction of about 2.7 x 
10~ 9 , including a phase space suppression factor of 0.75 in both cases. However, a 
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different choice of parameters could enhance the new physics contribution and make it 
competetive with the SM contribution. 

It is worth recalling that in the SM the decay matrix element is proportional 
to the square of a tree-level amplitude proportional to V c b ~ 4 x 1CT 2 arising from W 
exchange. As remarked in the Introduction, the best fit value for rf 2 is about twice as 
large as the SM prediction and it is challenging to understand this. For example, in the 
present model, second family Higgs exchange with mass Mh ~ 400 GeV and couplings 
in Eqs, 40|41f42|43[44 give a contribution to rf 2 which is completely negligible compared 



to the SM W exchange contribution. The corresponding charged Higgs exchange con- 
tributions are also expected to be suppressed compared to the SM. For example, the 
interaction DcbH^CRbi involves a coupling, 

Vcb ~ Vsb ~ e 2 d a d a h ~ 2.5 x 10" 3 , (45) 

which is again smaller than ~ 4 x 10~ 2 , with all charged Higgs couplings involving 
at least this suppression, and the charged Higgs mass being heavier than the W mass. 

Finally we remark that the model in [23] , as developed above, is based on A27 family 
symmetry combined with the EgSSM [32] which predicts three complete SUSY Higgs 
families as part of three complete 27 dimensional SUSY matter representations at the 
TeV scale (minus three right-handed neutrinos which get high see-saw scale masses since 
they carry no charges in this model). In addition there is a pair of SUSY doublets L, L 
in conjugate representations which form a TeV scale Dirac mass as required for GUT 
scale unification. These may be absent if the requirement of GUT scale unification is 
relaxed |33j. The E 6 SSM can also be tested via its prediction of a Z' N gauge boson with 
flavour conserving couplings 



4 Conclusion 

Following the anomalous like-sign dimuon charge asymmetry measured by the DO collab- 
oration at the Tevatron collider we have discussed the implications of large CP violation 
in B ds mixing for Supersymmetric (SUSY) Standard Models, focussing on those mod- 
els which allow a family symmetry and unification. For the Minimal Supersymmetric 
Standard Model (MSSM) we have seen from Eq{TT]that it is only possible to account for 
B s mixing and CP violation at the expense of large squark mixing which would require 
a new approach to family symmetry models. However, assuming such a framework, it 
seems very difficult to account for a significant amount of B d mixing and CP violation. 

In order to describe both B s and B d mixing and CP violation, as suggested by a 
recent global fit, we were led to consider SUSY models with Higgs fields transforming 
as triplets under a family symmetry. We have described a realistic such model based 
on A 27 family symmetry combined with the E 6 SSM in which tree-level exchange of the 
second Higgs family predicts B s and B d mixing and CP violation in the ratio h d /h s ~ 1/2 
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which is in good agreement with best fit point I in Eq|8j The model naturally suppresses 
flavour and CP violation in €k and the lepton sector, and is distinguished from other 
Higgs models by predicting Br(B s — > consistent with the SM prediction. 
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Appendix 

A The origin of phases 

In this Appendix we discuss the question of the origin of the phases in the flavour 



violating Yukawa couplings in Eq 36 Following the approach in [25] , we shall assume 
that CP is preserved in the high energy theory but is spontaneously broken by the flavon 
VEVs in Eqj27] whose phases can be restored as follows, 

(<fo> oc ( e^ 3 ) T , (<^) oc ( e^ 3+< ^ ) T , 

(023) oc (0 e iuJ2 e ^ 3+<fe ) ) T , (0i 23 ) oc ( e iwi e^+M e <(<"»+*») ) T , (46) 

where the phases Ui can be removed by SU(3) transformations but not in the A27 the- 
ory. Another difference between SU(3) and A27 is that the discrete symmetry allows 
nine distinct one dimensional representations [30], which, depending on the messenger 
representations, allows many more new operators than those given in Eqs.31, 32, corre- 
sponding to the different singlet contractions (3x3)i r where r = 1, ... 9. Here we restrict 
ourselves to A4 type messengers in the first three one dimensional representations which 
can be obtained from the products 3 x 3 as follows, 

1 = IT + 22 + 33, T = IT + w22 + w 2 33, l" = IT + w 2 22 + w33, (47) 

which are familiar from A4 [31] where u = exp(2ni/3). Allowing messengers in the 
1,1', 1" representations permits new operators corresponding to one messenger in each 
of the allowed one dimensional representations, where the invariant singlet is given by 
1 = V x 1". Thus the operators in Eqs |3~i|[32"1 need to be augmented by others of the more 
general form, (• ■ ■ )i(- • - )]/(• • • )i» appearing in all possible combinations. Assuming 
the messengers in the 1,1', 1" representations all have the same mass, the expansion 
parameters and predictions given previously will not change. In particular the flavour 
violating couplings will have their magnitudes unchanged from the values quoted in 
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Eqj36j but their phases will all be different from each other in a complicated way which 
depends on the order unity couplings which control the precise linear combinations of 
the different operators which contribute to these couplings. In the limit that only singlet 
operators are permitted (• • • )i(- • • )i(" " " )i it is easy to show that aig(yb s ) = aig(y s b) 



and axg(ybd) = arg(y^), even with the most general flavon VEVs in Eq,46, so the extra 
operators of the form (• • • )i(- • • )i'(- • • )i» are in fact necessary in order to allow new 
sources of CP violation in B s and mixing. 
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